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The stability of ferromagnetism at the surface at finite 
temperatures is investigated within the strongly correlated 
Hubbard model on a semi-infinite lattice. Due to the reduced 
surface coordination number the efi'ective Coulomb correla- 
tion is enhanced at the surface compared to the bulk. There- 
fore, within the well-known Stoner-picture of band ferromag- 
netism one would expect the magnetic stability at the surface 
to be enhanced as well. However, by taking electron corre- 
lations into account well beyond the Hartree-Fock (Stoner) 
level we find the opposite behavior: As a function of temper- 
ature the magnetization of the surface layer decreases faster 
than in the bulk. By varying the hopping integral within the 
surface layer this behavior becomes even more pronounced. 
A reduced hopping integral at the surface tends to destabilize 
surface ferromagnetism whereas the magnetic stability gets 
enhanced by an increased hopping integral. This behavior 
represents a pure correlation effect and can be understood in 
terms of general arguments which are based on exact results 
in the limit of strong Coulomb interaction. 



I. INTRODUCTION 



plies also for the magnetic transition metals Fe, Co, Ni, 
which are prototypical materials of the so-called band- 
magnetism of itinerant electrons. As it is well known, 
within a tight-binding scheme for the description of the 
bandstructure the reduced coordination number at the 
surface leads to a reduced effective bandwidth in the sur- 
face layer. As a consequence the effective correlation at 
the surface will be enhanced. Thus, on the basis of the 
well-known Stoner-picture of bandmagnetism one may 
intuitively expect the magnetic stability at the surface to 
be enhanced compared to the bulk. 

In this paper we investigate the magnetic stability at 
the surface in strongly correlated electron systems within 
the Hubbard model on a semi-infinite lattice. Based on 
the limit of strong Coulomb interaction we give general 
arguments concerning the magnetic stability at the sur- 
face. In addition we will show results of numerical evalu- 
ations within the spectral density approach to the semi- 
infinite Hubbard model. The spectral density approach 
is based on exact results about the general shape of the 
spectral density in the limit of strong Coulomb interac- 
tion. To emphasize the importance of correlation effects 
we will show, in addition, results obtained within the 
Hartree-Fock approximation (Stoner- model). 
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There is currently active interest in the influence of re- 
duced translational symmetry on the magnetic properties 
at surfaces, in ultra-thin films and in multi-layer struc- 
tures. One class of materials that is intensively stud- 
ied are the magnetic transition metals, within which the 
magnetically active electrons are itinerant. Theoretically 
it is possible to calculate the ground state properties like 
the magnetic moments at and near surfaces antLin thin 
films in great detail by use of ah initio methodad'0. How- 
ever, these approaches are strictly based on a Stoner-type 
picture of ferromagnetism and treat electron correlation 
effects that are responsible for the spontaneous magnetic 
order on a low level only. Further, an extension to finite 
temperatures states a difficult problem. 

To study the temperature behavior of spontaneous fer- 
romagnetism rather idealized model systems have proven 
to be a good starting point. What concerns the magnetic 
properties of systems with reduced translational sym- 
metry, the temperature behavior has been, up to now, 
almost exclusively studied withinrJocalized spin models 
(Ising model, Heisenberg model)0"El. In these systems the 
magnetic stability at the surface is found to be reduced 
compared to the bulk. The reduced magnetic stability re- 
sults simply from the reduced exchange coupling at the 
surface due to the lower coordination number. However, 
it is not clear at all to what extent this argument ap- 



II. THEORY 

To describe the geometry of the semi-infinite lattice 
each lattice vector is decomposed into two parts Ria = 
Ri + ■'"q- Throughout the paper Ri (Latin index) de- 
notes a lattice vector parallel to the film surface whereas 
r„ (Greek index) refers to the distance from the surface 
(a = 1: surface layer). Each lattice plane parallel to the 
surface contains N lattice sites. Within the layers we 
assume translational invariance. 

Using this notation the Hubbard model for the semi- 
infinite lattice reads: 
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Here Ciaa {^l^a) stands for the annihilation (creation) 
operator of an electron with spin a at the lattice site 



is the number operator. U denotes 



the on-site Coulomb matrix element and fj, the chemical 
potential. T^^ is the hopping integral between the lattice 
sites Ria and Rjp- 

Let us first discuss the kinetic part of the system 
{U — 0) that is treated within the tight-binding approxi- 
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mation. A two-dimensional Fourier transformation yields 
the corresponding dispersions 



k ij 



(2) 
(3) 



that can easily be calculated. Here and in the follow- 
ing k denotes a wave-vector from the underlying two- 
dimensional (surface) Brillouin zone, t is the hopping 
integral between nearest neighbor atoms. The on-site 
hopping integral Toq = T^"" gives the center of gravity of 
the dispersion of the a-th layer. The dispersions T^^ de- 
termine the Bloch density of states (BDOS) poa{E). The 
bandwidth of the BDOS is denoted by W. A moment 
analysis of the BDOS (aI"^ = / dEE'^paaiE)) shows 

the variance at the surface {a = 1) to be consider- 
ably reduced compared to the variance of the inner layers 

(2) 

Al (a = 2, . . . , oo). Mathematically, the square root 
of the variance {y^ a'^^) gives a measure of the effective 
bandwidth of the BDOS. In addition, for an asymmetric 
BDOS also the skewness Aq^'' is reduced at the surface 
compared to the bulk. For fcc-type semi-infinite lattices 
with surface orientations (100) and (111) that are con- 
sidered in the numerical evaluation (see below) we find: 

(Ai'VAf)('"°^ = 0.667, (A?Va1'^)("i) = 0.750, 



0.5, (Ai'VAp^)("i) = 0.625. 



(4) 



The corresponding surface coordination numbers are 
given by Zs^""^ = 8 and zi^^^-* = 9. These have to be 
compared to the bulk coordination number of the fee lat- 
tice zl^""^ = 12. 

The Hubbard model states a complex many-body 
problem, the rigorous solution of wMpb is possible up 
to now only for some limiting casesatll. In thin films 
and semi-infinite systems even more complications are 
introduced due to the reduced translational symmetry 
at the surface. The standard Hartree-Fock approxima- 
tion (HFA), which has been applied prepously for the 
investigation of thin film ferromagnetisinE3E3, is known 
to drastically overestimate the magnetic region in the 
phase diagram. This is clearly due to the fact that elec- 
tron correlations are treated on a low level within the 
HFA. Since fcrromagnetism in the Hubbard model is an 
effect of strong Coulomb interaction we require an ap- 
proximation scheme which takes correlation effects into 
account more reasonably. On the other hand, it must 
be simple enough to allow for a detailed study of phase 
transitions in semi-infinite systems. 

An approximation scheme which has proven to include 
the essential physics of spontaneous fcrromagnetism in 
the Hubbard model is the spectral density approach 
(SDA)yila. The SDA can be interpreted as a self- 
consistent extension of the W/U perturbation theory of 



Harris and Langeta in a straightforward way to inter- 
mediate coupling strength as well as finite temperatures. 
Despite its rather simple concept, the SDA leads to con- 
vincing results with regard to the qualitative behavior 
of the ferromagnetic solutions of the Hubbard model. A 
similar approach has been applied to a niultiband Hub- 
bard model with accurate results for the magnetic key 
quantities of the prototype ferromagnets Fe, Co, and 
NiO"ll3. Recently a generalization of the SDA has been 
proposed to deal with the mcwiifieations due to the re- 
duced translational symmetryn3~t3. In the following we 
give only a brief derivation of the SDA for a semi-infinite 
lattice. For fuKth i er . datails we refer the reader to previous 
publicationsym-ia. 

The basic quantity from which all relevant information 
on the system can be obtained is the retarded single- 
electron Green function 
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After a two-dimensional Fourier transformation one ob- 
tains from 

Gtf^iE) the spectral density S^j^{E) ^ 

— ilmG^^(i?), which represents the bare line-shape of a 
(direct, inverse) photoemission experiment. The diagonal 
elements of the Green function determine the spin- and 
layer-dependent quasiparticle density of states (QDOS): 
Paa{E) ~ —^lmGf^{E — n). Via an energy integration 
one immediately gets from Paa{E) the band occupations 
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(...) denotes the grand-canonical average and f-{E) is 
the Fermi function. Fcrromagnetism is indicated by a 
spin-asymmetry in the band occupations riaa leading to 
non-zero magnetizations ma — — Uai- The mean 
band occupation n and the mean magnetization m are 
given by n = i Y.aa "^aa and m = jY^a^a, which are, 
of course, identical to the respective bulk values (a 
oo). 

For the generalization of the SDA to systems with re- 
duced translational symmetry we adopt the local approx- 
imation for the self-energy {Tf^f^E) = Sij6°'f^T,aa{E)) 
which |-4ias been tested recently for a semi-infinite 
latticeca. The decisive step is to find a reasonable ansatz 
for the self-energy 'SacriE). Guided by the exactly solv- 
able atomic limit of vanishing hopping ii = 0, = 0) 
and by the findings of Harris and Langal3 in the strong- 
coupling limit (U/W ^ 1), a onedgole ansatz for the 
self-energy T,'^{E) can be motivatedEj'ES: 
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The spin- and layer-dependent parameters gf^, and 
g^^ are fixed by exploiting the equality between two al- 
ternative but exact representations for the moments of 
the Green function: 
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m- times 

Here [...,.. denotes the commutator (anticom- 

mutator). Eq. imposes rigorous sum rules on the 
Green function which have been recognized to state im- 
portant guideHnes when constructing approximate solu- 
tions for the Hubbard modelEJ. By comparing various 
approximation schemes it has been shown that the in- 
clusion of the first four sum rules is vital for a reason- 
able description of spontaneous ferromagnetisra_in the 
Hubbard model, especially at finite temperaturesEJ. Tak- 
ing the first four sum rules (m = 0, . . . , 3) into account 
to determine the free parameters of the SDA as 
for the self-energy yields a closed set of equationsE2r 
that has to be solved self-consistently. Note, that the 
self-energy incorporates the so-called bandshift Ba-a 
which is layer|-d ep j en cietit and consists of higher correla- 
tion functionalil't3E3"E3. Ba-a comes into play through 
the fourth sum rule (m = 3). Fortunately, the bandshift 
can pvartly hp pynrpgspH by the sclf-cnergy and the Green 
functionoEjcjCJ. A spin-asymmetry of the bandshift, 
which can occur within the self-consistent calculation in 
certain parameter regions, may generate or stabilize fer- 
romagnetic solutions in the strongly correlated Hubbard 
model. .For .strong Coulomb interaction the bandshift is 
foundE3EjO'E3 to be proportional to the kinetic energy of 

the system {na^a{\ - na-a)Ba~a ~ Ekin-a, U > W). 



III. RESULTS AND DISCUSSION 

The numerical evaluations have been done for a semi- 
infinite fee lattice. To study the influence of the reduced 
coordination number we have considered two different 
surface orientations (100) and (111). The mean band 
occupation is kept fixed at n = 1.4 and the Coulomb 
interaction is chosen to be J7 = iW which clearly refers 
to the strong coupling region. In order to avoid a mixture 
of charge transfer and electron correlation effects we have 
excluded charge transfer between the layers by requiring 
Ua = n, Va. This is achieved by adjusting the layer- 
dependent on-site hopping integral Tqq. 

Fig. 1^ shows the spontaneous magnetization of the 
semi-infinite fee lattice as a function of increasing tem- 
perature. At low temperatures the system shows full 
polarization (rn = 2 — n) and is uniformly magnetized 
{nia = m, Va) since we have explicitly excluded charge 
transfer between the layers. As the temperature increases 
the magnetization of the surface layer mi decreases faster 
than the bulk magnetization mi,. This behavior holds for 
both considered surface orientations, (100) and (111). 
The deviation from the bulk magnetization is almost 



completely restricted to the topmost layer. Note, that 
the influence of the surface is stronger for the more open 
(100) structure. Due to the coupling between the layers 
that is induced by the electron hopping a unique Curie 
temperature is found for the whole system. We have 
to conclude from Fig. |l], that within the strongly corre- 
lated Hubbard model the magnetic order is less stable 
at the surface than in the bulk, although the effective 
correlation is enhanced at the surface. This observation 
clearly contradicts the expectation on the basis of the 
well known Stoner-picture of band ferromagnetism. In- 
deed, repeating the same calculations within the HFA we 
find the magnetic order to be more stable at the surface 
compared to the bulk. For very strong Coulomb inter- 
action the layer magnetizations within the HFA become 
almost independent of the layer index. 

The reduced magnetic stability at the surface in 
strongly correlated electron systems can be understood in 
terms of general arguments by use of the moment analy- 
sis of the BDOS of the semi- infinite lattice (see Eq. (^). 
The arguments are based on the magnetic behavior of 
bulk systems with full translatiouaJ symmetry, which 
have been studied in very detailExES. Both the reduced 



variance < A^^'') as well as the reduced skewness 
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(aI^'' < A^''') at the surface lead to the same trend: 



,(3) 
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FIG. 1. Magnetization nia for a semi-infinite fee lattiee as 
a function of the reduced temperature T/Tc- a = 1 cor- 
responds to the surface layer, (a) (100) surface; (b) (111) 
surface. Further parameters: U/W = 3; n = 1.4 



(i) Within various approximation . schemesr"H""H"i as 
well as DMFT-QMC calculationsBGa it has been estab- 
lished, that a large skewness of the BDOS tends to stabi- 
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lize spontaneous ferromagnetism in the Hubbard model. 
Since the skewness of the surface layer BDOS is reduced 
compared to the inner layers this explains the trend to a 
reduced magnetic stability at the surface. 

(ii) In order to understand the influence of the reduced 
variance of the BDOS at the surface on the magnetic sta- 
bility we have plotted in Fig. ^jthe Curie temperature Tc 
for an fee bulk system in two different ways. Fig. ||(a) 
shows kBTc/W as a function of the effective correlation 
U/W. This corresponds to the situation where the band- 
width W of the uncorrelated Bloch band is kept fixed 
while the Coulomb interaction U is varied. Within the 
SDA the Curie temperature steeply increases as soon as 
U/W exceeds a critical lower bound. But already for 
U/W sa 1 — 2 the Curie temperature starts to satu- 
rate and approaches a finite value as U — >• cx). This 
behavior can be understood by inspecting the quasiparti- 
cle spectrum. In the strongly correlated Hubbard model 
the quasiparticle spectrum splits into a high and a low 
energy subband (" HubbauL-bands" ; for a detailed dis- 
cussion see, e. g., Refs.BEl). Both subbands are sepa- 
rated by an energy amount of approximately U. In the 
ferromagnetic phase there is an additional spin splitting 
within each subband between the majority and the mi- 
nority spin direction. According to the results of Har- 
ris and LangeliS, the spin splitting is - for large U - 
determined by the spin-asymmetry in the kingtic .en.- 
ergy Eki„ = Ekin.i + Ekirui of t^^e systemEjMH'B. 
This can be seen, e. g., for the spin-splitting A2.ex be- 
tween the centers of gravity (T'2ct) of the upper Hubbard- 
bands within which the chemical potential fi is located 
for t|bepa»0|ce, than half-filled case (n > 1). Here one 



findsE 
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(9) 



The kinetic energy, however, becomes independent of U 
in the limit of strong Coulomb interaction and the band- 
width W remains the only relevant energy scale with re- 
spect to the magnetic properties of the system. There- 
fore, the spin splitting (|^) as well as the Curie tempera- 
ture saturate as a function of increasing U (Fig. ||(a)). 
The same behavior has been observed within . soyo ral 
strong coupling approaches to the^Hubbard modelEacJ as 
wefl as DMFT-QMC calculations^. Within the HFA the 
spectrum shows a completely different behavior. Here, 
the Hubbard splitting is absent and the spin splitting is 
given by Um. Both, the spin splitting and the Curie tem- 
perature increase linearly with U in the strong Coupling 
limit as can be seen in Fig. |2|(a). In Fig. ||(b) we have 
plotted ksTc/U as a function of the effective correlation 
U /W . This means that the Bloch bandwidth W is varied 
while the Coulomb interaction U is kept fixed. Note, that 
this corresponds to the situation in semi-infinite lattices 
as well as in thin films where U is uniform for all layers. 



the effective bandwidth ^/a^) of the BDOS being, 
however, different at the surface and in the bulk. As we 
have discussed above, in the strong coupling regime the 
energy scale being important for the magnetic properties 
is given by W . Therefore, for fixed U the Curie tem- 
perature within the SDA decreases proportional to 1/W 
as U/W 00 (see Fig. |(b)). In the limit ly = we 
find a vanishing Curie temperature, which reproduces the 
exact result in the so-called "atomic limit" of vanishing 
hopping. Again, the Hartree-Fock solution fails to de- 
scribe the strong coupling behavior correctly, giving a fi- 
nite Curie temperature even for W — Q. Since the energy 
scale determining the magnetic properties is given by U 
within the HFA the Curie temperature saturates at a fi- 
nite value as ly — > in Fig. ^(b). By applying the results 
of Fig. |(b) to the semi-infinite system it is clear that in 
the limit of strong Coulomb interaction the reduced ef- 
fective bandwidth in the surface layer (y^A^ < A<^') 
enhances the effective correlation but reduces the mag- 
netic stability at the surface. In addition, Fig. ^(b) ex- 
plains the qualitatively different behavior concerning the 
magnetic stability at the surface within the HFA and the 
SDA. 
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U/W 

FIG. 2. Curie temperature Tc as a function of the effec- 
tive correlation U/W for an fee lattice (bulk calculation), (a) 
constant bandwidth W of the BDOS . (b) constant Coulomb 
interaction U. fcs is the Boltzmann constant. Further pa- 
rameters: n = 1.4 

We want to stress that both arguments (i) and (ii) that 
explain the reduced magnetic stability at the surface in 
strongly correlated electron systems are not restricted or 
linked to the SDA. Qualitatively, the same behavior is to 
be expected within all theories that correctly reproduce 
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the exact results of the W/U perturbation theory of Har- 
ris and LangetJ. It is clear that argument (ii) holds only 
for strong Coulomb interaction. In the case of moderate 
coupling strengths the effective correlation U/W has to 
exceed a lower critical bound for spontaneous ferromag- 
netism to occur. 




FIG. 3. ^(^) (see Eq. Q) calculated within the HFA 
and the SDA as a function of the effective correlation U/W for 
an fee lattice (bulk calculation). Further parameter: n = 1.4. 

Let us discuss the boundary between these two regimes 
(moderate coupling <-> strong coupling behavior) in more 
detail. Due to the scaling pmperty of the Hubbard model 
the following relation holdsLj: 



Tc{U,rW)^rTc{-,W), 
r 



r > 



(10) 



Eq. |l^ expresses the fact that for constant U/W the en- 
ergy scale of the Hubbard model is determined by W . A 
Taylor expansion of the right hand side of (|l0|) at r = 1 
up to the first order yields: 



Tc{U,rW)^Tc{U,W) 



U fdTc 



(1-r) 



(11) 



For -M-i^S^) > 1 a reduction of the bandwidth W 
by a factor r will increase the Curie temperature, for 
^(^^) < 1 the Curie temperature will be reduced. 

In Fig. § numerical results for ^(^^) are plotted as 
a function of U/W for both the SDA and the HFA. 
^(^^), that is defined in the ferromagnetic phase only, 
diverges as U/W approaches the lower critical corre- 
lation. Thus, for moderate coupling strength we find 
^(^) > 1 within the HFA and the SDA. However, 

within the SDA ^(^^) becomes smaller than 1 very 
rapidly with increasing U/W and vanishes as U/W oo. 
Thus, by inspecting ^(^^) we can distinguish between 
two different regimes: For moderate correlation we ex- 
pect the magnetic stability at the surface to be enhanced, 



whereas it will be reduced for strong Coulomb interac- 
tion. Again, a similar behavior is to be expected within 
other strong coupling approaches to the Hubbard model 
as well. 




FIG. 4. Bulk magnetization nib and surface layer magne- 
tization mi as a function of the reduced temperature T /Tc 
for a semi-infinite fee lattice with a (111) surface. The hop- 
ping integral within the first layer has been varied according 
to Tfe" = Toi + rf7i|(fe), e > 0. (a) e = 1.0, 0.9, 0.8, 0.7. 
(b) e = 1.0, 1.1, 1.2, 1.3. The thick lines correspond to the 
uniform situation e = 1.0. Further parameters: U/W = 3; 
n = 1.4 

To further investigate the magnetic stability at the sur- 
face we have modified the hopping integral within the 
surface layer of the semi- infinite fcc(lll) lattice accord- 
ing to Tl^ = Toi + et7||(fc), e > 0. In Fig. | the corre- 
sponding layer magnetizations are shown as a function of 
the reduced temperature T/T^""^. We have considered 
deviations from the unperturbed hopping e = 1.0 within 
30% only. A reduced hopping integral in the surface layer 
(e < 1, Fig. ^(a)) results in a reduced magnetic stability 
at the surface. As in the uniform situation (e = 1, see also 
Fig. l|) we find a unique Curie temperature Tc — rp"""' 
for all e < 1 due to the coupling between the surface layer 
and the bulk. As the hopping at the surface is enhanced 
compared to the bulk (e > 1, Fig. ^(b)), the magnetic sta- 
bility at the surface increases. Note that already a 30% 
enhancement of the surface hopping integral (e — 1.3) 
leads to an enhanced surface Curie temperature T^^^ 



which lies clearly above T^""^. 



For r^""= 



<T <Tc 



the polarization of the surface layer induces a finite po- 
larization in the adjacent layers. However, the induced 
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magnetization decreases very fast with increasing dis- 
tance from the surface. In this sense we find two sep- 
arate Curie temperatures in the system. The bulk Curie 
temperature is, of course, not affected by the enhanced 
surface Curie temperature. Quahtatively, the influence 
of the modified hopping at the surface on the magnetic 
stability can be understood by argument (ii) that was in- 
troduced above. A reduced (enhanced) hopping integral 
at the surface leads to a reduction (enhancement) of the 
effective bandwidth of the surface layer BDOS and, there- 
fore, to a reduced (enhanced) magnetic stability. Again, 
the opposite behavior is found within the HFA. 

IV. CONCLUSION 

In conclusion we have investigated the magnetic sta- 
bility at the surface in strongly correlated electron sys- 
tems within the semi-infinite Hubbard model. For strong 
Coulomb correlation we find the magnetic stability at the 
surface to be reduced compared to the bulk. This behav- 
ior represents a pure correlation effect that cannot be 
explained within the well-known Stoner-picture of band- 
magnetism which is justified for weak Coulomb interac- 
tion only. However, based on exact results concerning the 
general shape of the spectral density in the strong cou- 
pling limit the reduced magnetic stability at the surface 
can be understood by a moment analysis of the BDOS. 
It has been shown that both the reduced variance (effec- 
tive bandwidth) as well as the reduced skewness of the 
surface layer BDOS tend to reduce the magnetic stability 
at the surface. 

By modifying the nearest neighbor hopping integral 
within the surface layer we have investigated the mag- 
netic stability at the surface in more detail. Again, the 
same trend is observed. A reduced (enhanced) hopping 
integral at the surface leads to a reduced (enhanced) mag- 
netic stability at the surface. For an fcc(lll) geometry 
already a 30% enhancement of the surface hopping in- 
tegral results in a surface Curie temperature which lies 
well above the bulk Curie temperature. 
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